DESY 08-161 
HU-EP-08/55 
SFB/CPP-08-89 



A comparison of the cut-off effects for 
Twisted Mass, Overlap and Creutz fermions 
at tree-level of Perturbation Theory 

KRZYSZTOF CICHY 
Adam Mickiewicz University, Faculty of Physics, 
Umultowska 85, 6I-6I4 Poznan, Poland 

JENIFER GONZALEZ LOPEZ 

Humboldt-Universitdt zu Berlin, Institut fiir Physik, 
Newtonstrasse 15, 12489 Berlin, Germany 
and 

DESY Zeuthen, Platanenallee 6, D- 15738 Zeuthen, Germany 

AGNIESZKA KUJAWA 
Adam Mickiewicz University, Faculty of Physics, 
Umultowska 85, 6I-6I4 Poznan, Poland 

Abstract 

In this paper we investigate the cutoff effects at tree-level of pertur- 
bation theory for three different lattice regularizations of fermions - 
maximally twisted mass Wilson, overlap and Creutz fermions. We 
show that all three kinds of fermions exhibit the expected O^a"^) scal- 
ing behaviour in the lattice spacing. Moreover, the size of these cutoff 
effects for the considered quantities i.e. the pseudoscalar correlation 
function Cps, the mass mps and the decay constant fps is compara- 
ble for all of them. 
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1 Introduction 



The main goal of Lattice Field Theory is to study the non-perturbative as- 
pects of quantum field theories, in particular Quantum Chromodynamics 
(QCD). For example, Lattice QCD is a regularization of QCD which con- 
sists in putting the theory on a four dimensional lattice (discretization) with 
lattice spacing a, whose inverse is the ultraviolet cutoff of the theory. The 
discretization of bosons is relatively straightforward, but when one tries to 
discretize fermions in the naive way, the notorious fermion doubling problem 
emerges - instead of one fermion in the continuum limit, one has as many as 
2^ fermions, where d is the space-time dimensionality. As it was originally 
proposed by Wilson [IJ, the doubling problem can be solved if a different 
fermionic discretization is chosen, the so called Wilson fermions. After Wil- 
son's proposal, many alternative fermion discretizations which remove the 
doubling problem have been suggested and are still appearing. However, 
as stated by the Nielsen-Ninomiya theorem [2], new problems will always 
appear when removing the doublers; in order to eliminate the fermion dou- 
bling problem one has to pay the price of either explicitely breaking chiral 
symmetry (even in the massless limit), or giving up locality or translational 
invariance. 

Much of the effort of lattice QCD goes into finding a lattice theory without 
doublers which keeps the largest possible number of symmetries, and at the 
same time reaches the continuum limit as fast as possible (the dependence 
on the inverse cutoff, a, is as small as possible e.g. 0{a'^) leading cutoff 
dependence is better than 0{a)). 

In this paper, we investigate the cutoff effects at tree-level of perturbation 
theory of three different discretizations of fermions - twisted mass Wilson 
fermions at maximal twist (MTM), overlap fermions and Creutz fermions, 
at a fixed value of the physical quark mass. We have presented a similar 
analysis for a different value of the quark mass in [3j. Here we compare the 
results. 

The MTM fermions [1], [H] are relatively cheap to simulate and they are 
by now a widely used fermion discretization. Although similar to Wilson 
fermions, they retain a subgroup of chiral symmetry which guarantees au- 
tomatic 0{a) improvement, i.e. 0{a?) leading cutoff effects. The price to 
pay to have a residual chiral symmetry, is to break a subgroup of the isospin 
symmetry transformatior0. 

It has been shown by Ginsparg and Wilson [7] that there is a way to 
preserve chiral symmetry on the lattice, even without the doubler modes, if 

"'^For a recent review of twisted mass fermions see [1]. 
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the corresponding Dirac operator obeys a relation now called the Ginsparg- 
Wilson relation. It is a non-standard realization of chiral symmetry [S], be- 
cause the Dirac operator no longer anticommutes with 75 at non-zero lattice 
spacing, but it only anticommutes with a lattice modified version of 75. A 
particularly simple form of a Dirac operator that obeys the Ginsparg- Wilson 
relation has been found by Neuberger [9j. The main disadvantage of overlap 
fermions is that they are much more costly to simulate - by a factor of 30-120 
in comparison with MTM fermion^. 

The recently proposed Creutz fermions [T2] represent the class of minimally- 
doubled fermionqj. They describe two flavours of quarks and preserve exact 
chiral symmetry. However, they break a number of discrete symmetries and 
isospin symmetry and this would make their simulation very difficultQ. 

2 Setup 

2.1 Correlation functions 

In order to investigate the cutoff dependence of the pseudoscalar meson mass 
and decay constant we have to first calculate the correlation function cor- 
responding to this mesorif]. Despite the fact that we work at tree-level of 
perturbation theory, we will refer to the pseudoscalar meson as the 'pion'. 
The charged pions are desribed by the following interpolating operator: 

r^{x) =V\x)TiV\x), (1) 

where the pseudoscalar density P"(a;) = ^/'(x)75^^/'(x) (for a = 1,2,3) and 
r" are the Pauli matrices. 

The time dependence of the correlation function Cps{t) is thus given by: 

Cps{t) = -Y.{0\V^i^)V-m0). (2) 

X 

Performing all the possible Wick's contractions, one obtains the dependence 
of the pion correlation function on the quark propagator, S'^(p), given by 

Cpsit) = ^ E E E E e'^'^-''^^'s,{p,p,)s;{p,p',). (3) 

^For a review of overlap fermions see e.g. [10 . For a comparison with twisted mass 
fermions see [TT| . 

Other examples of minimally-doubled fermions were given by Karsten [13] and Wilczek 

m- 

^For a discussion of this aspect of Creutz fermions see [T3] . 

^For a pedagogical introduction to the methods we have used in this work see [IB] . 
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Nc denotes the number of colours and the number of Dirac components. 
L = aN is the physical extent of the lattice in the spatial directions {N is 
the number of lattice sites in all spatial directions) and T = aN^ the physical 
extent in the temporal direction (A^4 - the number of lattice points in the 
time direction). The possible choices of the index fi will be explained below 
(see Eq. f ill I) ). The numerical computation of correlation functions consists 
in directly evaluating the expression ([3]). 



2.2 Quark propagators 

We present in this section the analytical expressions of the quark propaga- 
tors, in momentum space and at tree-level of perturbation theory, for the 
three kinds of lattice fermions considered in our analysis. 



Wilson twisted mass fermions 



where: 

p^, = - sin(ap^), p^, = - sin(^), M{p) = ruo + ^ J2 P% (5) 



1 and 1 J are the identity matrices in Dirac and flavour space, respectively, [iq 
is the twisted quark mass and mo the untwisted quark mass. The maximal 
twist setup consists in setting the untwisted mass to zerd^ such that the 
quark mass is only given by the twisted mass. 



Overlap fermions 



-z(l-fi)F(p)-V2^^7^ + _A4(p)]l 



(1 - ^ypip)-^ E.Pl + Mipy ' 



where: 

F{p) = 1 + ^T.pIpI, (7) 



2 IJ,<U 



m is the bare overlap quark mass. 

^This can be done exactly only at tree-level. A fine tuning is required in the interacting 
theory. 
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Creutz fermions 



Scip) 



where: 



(9) 
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7 = a" 7, mo - bare quark mass, C - lattice geometry parameter, R - C- 
dependent normalization factor needed to obtain the correct continuum limit. 
We consider two values of C: C = 3/^10 {R = 2)0 and C = 3/VU {R = 

We also consider a modification of Creutz's action suggested by Borici 
|17j . We call the corresponding fermions the 'Borici fermions' and the quark 
propagator for them is: 



^b(p) 



-i G^{ap)jf, + mo 1 



Eu G^(ap)2 + mg 



(10) 



where the functions Gfj,{ap) are: 



Gi{ap) 
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Quark propagator decomposition. All of the quark propagators are 
matrix expressions that can be decomposed in terms of the gamma matrices 
and the identity matrix: 



sip) = Suip)i + J2s,iph,, 



'IV. 



where /i = 1,2, 3,4 for overlap and Creutz fermions and /i = 1,2,3,4,5 for 
twisted mass fermions. 



^This value corresponds to the hypercubic lattice. 

^This value corresponds to a highly symmetric lattice geometry, which is the 4- 
dimensional analogue of graphene structure. 
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N'^Cpsit/N = 4) 


a 


b 


c 


MTM 


0.00996934 


0.00170143 


0.00002268 


OVERLAP 


0.00996934 


-0.00021268 


-0.00006924 


BORICI 


0.00996934 


0.00329653 


-0.00116956 


CREUTZ - C = 3/VW 


0.00996934 


0.00499799 


0.000201048 


CREUTZ - C = 3/714 


0.00996934 


0.00412066 


-0.00143986 



Table 1: Fit coefficients for the pseudoscalar correlation function. 



3 Scaling tests 

In this section we present the scaling tests performed on the pseudoscalar 
correlation functions, masses and decay constants. We employ the following 
strategy; we fix Nm = 0.8 (where m is the bare quark mass in lattice units) 
and calculate the correlators towards the continuum limit (a — > 0). At tree 
level of perturbation theory this is equivalent to the limit N oo. The time 
extent is always set to be larger than and an integer multiple of the spatial 
extent. 

In Fig. [T]we show the correlation function at a fixed physical time t/N = 
4, which is large enough to allow for a reliable extraction of the ground 
state contribution. To compare different fermion discretizations, we extract 
the coefficients. Table [H of the fitting curves shown in Fig. [H We use the 
following form of the fitting function: 

N'Cps = a + b^^+c^^. (12) 

Fig. [2] shows the pseudoscalar mass and in Table [2] we have gathered the 
coefficients of the following fit: 

Nmps = a + b^ + c^^. (13) 

Fig. [3] presents the pseudoscalar decay constant and Table [3] the coeffi- 
cients of the following ffi: 

iV/p5 = a + &^ + c^. (14) 

All types of fermions show the expected behaviour in the lattice spacing - 
O(a^) scaling violations. This is due to the exact chiral symmetry for overlap 
and Creutz fermions, and to the residual chiral symmetry for MTM. 

The continuum limit for each observable is always the same for every 
discretization (and the expected one for the mass at tree-level of perturba- 
tion theory), thus providing a ffist check of consistency of the corresponding 
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Figure 1: Cutoff effects and continuum limit of the pseudoscalar correlation 
function. 
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Figure 2: Cutoff effects and continuum limit of the pseudoscalar mass. 
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Figure 3: Cutoff effects and continuum limit of the pion decay constant. 



Nmps 


a 


b 


c 


MTM 


1.6 


-0.042667 


0.003045 


OVERLAP 


1.6 


0.085333 


0.008182 


BORICI 


1.6 


-0.202667 


0.058999 


CREUTZ - C = 3/VTO 


1.6 


-0.032000 


-0.106063 


CREUTZ - C = 3/VT4 


1.6 


-0.200000 


0.029668 



Table 2: Fit coefficients for the pseudoscalar mass. 



Nfps 


a 


b 


c 


MTM 


2.73861 


0.109545 


-0.004307 


OVERLAP 


2.73861 


0.219089 


0.028606 


BORICI 


2.73861 


-0.136931 


-0.027123 


CREUTZ -C = 3/VlO 


2.73861 


0.593370 


-0.388244 


CREUTZ - C = 3/\/l4 


2.73861 


-0.015977 


-0.195636 



Table 3: Fit coefficients for the pseudoscalar decay constant. 
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lattice regularizations here analyzed. The magnitude of the O(a^) effects 
is, however, very different for different discretizations and depends on the 
observable under consideration. For example, for the correlation function at 
a fixed physical time, the smallest effects are exhibited by overlap fermions 
and the largest by Creutz fermions with C = 3/a/TO. For the pion mass, 
however, the O(a^) scaling violations are the smallest for Creutz fermions 
with C = 3/vTO and the largest for Borici fermions. Therefore, there are 
no definite conclusions, from this scaling test at tree-level, of which fermions 
exhibit the smallest O(a^) effects. The only clear regularity that we ob- 
serve is that the disretization errors for twisted mass fermions at maximal 
twist are rather small for all observables that we have considered. Even the 
recently proposed Creutz fermions and their modification by Borici, which 
break a number of important discrete symmetries, do not suffer from very 
large O(a^) scaling violations at tree- level of perturbation theory and thus 
can not be excluded from this point of view. 

We have compared the results for the scaling behaviour here presented 
with the ones discussed in [3], where the same study for a different value of 
the quark mass {Nm = 0.5) was performed. We observe that, as expected, 
for all the quantities which have a well-defined continuum limit, the relative 
discretization errors (ratio of the coefficient of the O(a^) effects with respect 
to the continuum value) decrease when decreasing the quark mass indepen- 
dently of the action considered, while the difference of the relative discretiza- 
tion errors can vary between the actions considered here when changing the 
quark mass. 

4 Conclusion 

We have performed a scaling test of three different lattice fermion regulariza- 
tions at tree-level of perturbation theory; the widely used twisted mass and 
overlap fermions and also the recently proposed minimally-doubled Creutz 
fermions. All these discretizations lead to the same continuum limit and are 
C(a)-improved, but the relative sizes of the C(a^) effects depend strongly 
on the observable we choose for the analysis. Therefore, we can not exclude 
or put a preference, from the tree-level study of the lattice artifacts, on any 
particular fermion discretization of the three here considered. It will there- 
fore be interesting to test these discretizations in the interacting theory in 
practical simulations. 
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